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ABSTRACT 

^  Kinematic  precision  is  affected  by  errors  which  are  the  result  of  either 
intentional  adjustments  or  accidental  defects  in  manufacturing  and  assembly  of 
gear  trains.  This  paper  explains  a  general  method  for  the  determination  of 
kinematic  precision  of  gear  trains.  The  general  method  is  based  on  the  exact 
kinematic  relations  for  the  contact  point  motions  of  the  gear  tooth  surfaces 
under  the  influence  of  errors.  An  approximate  method  is  also  explained. 

Example  applications  of  the  general  and  approximate  methods  are  demon¬ 
strated  for  gear  trains  consisting  of  involute  (spur  and  helical)  gears, 
circular-arc  (Wi ldhaber-Novikov)  gears,  and  spiral-bevel  gears.  Gear  noise 
measurements  from  a  helicopter  transmission  are  presented  and  discussed  with 
relation  to  the  kinematic  precision  theory. 
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INTRODUCTION 


Transmission  error  is  a  measure  of  the  kinematic  precision  of  gear  trains. 
Transmission  error  is  defined  as  the  departure  of  the  meshed  gear  pair  (or  en¬ 
tire  gear  train)  from  a  constant  ratio  of  angular  motions  as  defined  by  the 
ratio  of  tooth  numbers.  It  is  true  that  in  a  gear  pair  each  gear  has  a  whole 
number  of  teeth  and  this  defines  the  nominal  ratio  of  angular  positions  between 
the  two.  But  the  instantaneous  ratio  during  a  meshing  cycle  can  vary  slightly 
from  the  nominal  ratio.  Transmission  error  is  the  measure  of  instantaneous 
variation  from  the  ideal  nominal  value.  The  precision  of  gears  was  investi¬ 
gated  by  Litvin  [1]*,  Litvin  and  Gutman  [2],  and  Michalec  [3], 

When  the  mating  teeth  in  a. gear  train  have  profiles  that  transmit  motion 
having  no  error,  they  are  said  to  be  conjugate  pairs.  In  theory,  it  is  pos¬ 
sible  to  select  an  arbitrary  shape  for  a  driving  tooth  and  then  to  find  a  pro¬ 
file  for  the  driven  tooth  which  will  give  conjugate  action.  Involute  gear 
teeth  happen  to  have  the  same  form  for  driving  and  driven  member  teeth.  A 
benefit  of  the  involute  form  is  that  small  errors  in  center  distance  between 
gears  will  not  produce  transmission  errors.  This  is  not  true  for  other  tooth 
profile  forms  [4,  5]. 

In  general,  there  are  many  causes  for  transmission  error,  and  they  cannot 
be  avoided  in  practice.  Such  things  as  shaft  misalignment,  profile  error, 
tooth  deflections  under  load,  mounting  location  errors,  and  gear  support 
deflections  may  combine  to  cause  transmission  error.  The  effects  of  trans¬ 
mission  error  are  most  often  harmful.  These  are  high  vibration  and  noise, 
pitting  and  scoring  of  gear  teeth,  and  reduced  reliability  of  the  gear  train. 
Sometimes  it  is  beneficial  to  introduce  small  intentional  errors  into  the  gear 


♦Numbers  in  square  brackets  denote  references. 


tooth  profile  to  compensate  for  the  probable  occurrence  of  accidental  and  una¬ 
voidable  errors  in  the  assembled  and  operating  gear  train.  Tip  relief  to  re¬ 
duce  dynamic  loading  and  combined  mismatch  in  spiral-bevel  gears  to  reduce  mis¬ 
alignment  sensitivity  are  two  examples  of  intentional  errors  which  are  bene¬ 
ficial. 

Baxter  has  studied  the  effect  of  various  types  of  misalignment  on  tooth 
contact  in  bevel  and  hypoid  gears  [6],  Townsend,  Coy,  and  Hatvani  have 
examined  gear  train  noise  as  a  test  of  its  precision  during  an  intentional 
loss  of  lubricant  destruction  test  [7]. 

Of  course,  the  effect  of  all  errors  (intentional  or  otherwise)  on 
transmission  error  must  be  predicted  analytically  if  the  gear  design  process 
is  to  remain  rational  and  not  collapse  into  a  confusing  heap  of  empiricism. 

Errors  of  manufacturing  and  assemblage  of  gears  induce  kinematic  errors 
in  gear-drives  which  may  be  presented  by  the  following  function: 

(1) 

Here  is  the  angle  of  rotation  of  the  driving  gear,  1, 

aQ  =  Uqj,  aq2,  .  .  .)  (2) 

is  the  vector  of  errors,  and 

2  ■  *2  “  *2  (3) 

is  the  kinematic  error  of  the  gear  drive,  represented  as  the  difference  between  the 
theoretical  and  actual  angles  of  rotation  of  driven  gear  2. 

In  this  paper  two  methods  to  determine  Function  (1)  are  presented:  (a)  a 
numerical  method  for  computer  solution  and  (b)  an  approximate  geometric  solution 
which  leads  to  simple,  accurate  results  in  an  analytical  form. 


NOMENCLATURE 


shortest  distance  between  gear  axes  of  rotation 
eccentricity  vector  of  gear  i 
axial  settings  of  gears  (see  Fig.  2) 
angular  velocity  ratio,  gear  i,  gear  j 
number  of  teeth  on  gear  i 

vector  function  representing  unit  normals  of  surface  of  gear  i 
in  fixed  coordinate  system  (f) 

change  in  unit  normal  vector  due  to  errors  in  gear  i 

change  in  unit  normal  vector  due  to  point  motion  relative  to 
gear  i 

change  in  unit  normal  vector  due  to  transfer  point  motion  with 
gear  i 

absolute  velocity  of  tip  of  unit  normal  vector  of  surface  i 

similar  to  dn^  and  dn^  but  velocities  rather  than 
_r  _tr 

displacements 

center  of  rotation  of  gear  i 

geometrical  center  of  gear  i 

pitch  point 

vector  of  errors 

components  of  vector  of  errors 

vector  from  origin  to  axis  of  gear  rotation 

base  circle  radius  of  gear  i 

vector  function  representing  surface  of  gear  i  in  fixed 
coordinate  system  (f) 


velocity  of  contact  point 
coordinate  system  i 

coordinate  system  rigidly  connected  to  frame 
absolute  displacement  of  contact  point  of  gear  i 

displacement  of  contact  point  due  to  errors  in  gear  i 

displacement  of  contact  point  relative  to  gear  i 

displacement  of  contact  point  in  transfer  motion  with  gear  i 
tangent  plane 

surface  coordinate  of  gear  i  surface 
similar  to  ds^gi  dstr^’  but  vel°c'it''es 

rather  than  displacements 

angular  position  of  eccentricity  vector  of  gear  i 

initial  angular  position  of  eccentricity  vector  of  gear  i 

rotation  vector  representing  position  change  of  gear  axis  of 
rotation 

surface  coordinate  of  gear  surface  i 

change  in  kinematic  error  function  as  measured  on  shaft  of 
gear  i 

proportionality  factor 
radius  vector 
surface  i 

angle  of  rotation  of  gear  1 
actual  angle  of  rotation  of  gear  2 

theoretical  angle  of  rotation  of  gear  2 


kinematic  error  function 


L. 

Tq  pressure  angle 

angular  velocity  of  gear  i 


THEORY  AND  EXACT  SOLUTION  METHOD  FOR  KINEMATIC  PRECISION 
In  the  process  of  motion  the  tooth  surfaces  of  two  gears,  and  £ 2 
(fig.  1),  are  in  tangency  if  the  following  equations  are  satisfied: 


.(D 


rf^(u2,e2’ 


,(D 


(2) 


( ^2  »®2  2  ^ 


Here  r|^  is  the  position  vector  of  the  contact  point  on  gear  i;  is 

the  surface  unit  normal  vector  at  the  contact  point  M;  u^,  and  are 
the  surface  coordinates  of  the  gear  surfaces;  and  ^  is  the  angle  of 
rotation  of  gear  i.  Subscript  f  denotes  a  coordinate  system  which  is  rigidly 
connected  to  the  frame. 

For  a  gearset  with  kinematic  errors,  represented  Dy  aQ^  and  aQ^,  con¬ 
ditions  for  tangency  may  be  expressed  as 


ff  ( 0^ ,©2  ,<p 2 > aQj )  =  r^  ^  ^2*^2 ’  ^2  ^ 
n|  ^  yf  j » aQj  )  =  nj.  ^  ( u^ ,©2  ,<p2 , aQ^ ) 


(< 
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Equations  (6)  and  (7)  yield  the  functions 

92(<?1»AQ1,aQ2)  =  v  2^  1 )  +  av2(»i»aQi»aQ2) 
Ui i *aQi , aQo )  ;  »aQi *aQo) 


(< 


The  functions 


r^1,(ui>«i)  ;  u.(v 1,aQ1,aQ2),  •1(«1,aQ1,aQ2)  (i  =  1,2)  (10) 

represent  the  path  of  the  contact  point  on  gear  surface  t.  corresponding  to  the 
meshing  of  gears  with  errors  of  manufacturing  and  assembly.  Functions 

rfi)(u?(*l)»ei(*l))  ;  ui(vlK  ei<*l>  UD 

represent  the  path  of  the  contact  point  on  gear  surface  i^  corresponding  to 
meshing  without  errors.  Comparison  of  functions  (10)  and  (11)  yields  the  change 
of  the  contact  point  path  induced  by  errors. 

Consider  the  solution  of  equations  (4)  and  (5)  and  (6)  and  (7).  Vector 
equations  (4)  and  (5)  yield  only  five  independent  scalar  equations  since 
n(l)|  =  |n(2)|  =  !•  These  equations  may  be  presented  as 

f  j  ( j  tV  j ,<p^)  =0  (j  —  1,2,.. .,5)  (12) 

It  is  assumed  that 

|fl>f2»f3,f4,f5}  e  d 

The  symbol  d  denotes  that  functions  f .  have  continuous  partial  deriva- 

v 

tives  of  the  first  order  (at  least)  by  all  arguments. 

It  is  assumed  that  equation  system  (12)  is  satisfied  by  a  set  of  para¬ 
meters 


p(l) 


,9 


(1) 

2 


) 


(13) 


and  that  surfaces  i^  and  ip  are  in  tangency  at  point  M  .  Surfaces  i^ 
and  will  be  in  point  contact  in  the  neighborhood  of  Mq  if  by  the  set  of 
parameters  the  following  Jacobian  differs  from  zero: 


l»t 2»^3»^4»^5) 
D(Uj ,©^ » ^2  *^2  2  ^ 


3fl 

3  f  1 

3  f  1 

3f  1 

3  f  1 

3U^ 

3©^ 

3U£ 

3«2 

S<?2 

3f5 

3;5 

3^S 

af5 

3f5 

3  U^ 

3©^ 

3U2 

3©2 

3<p? 

*  0 


(14) 


If  inequality  (14)  is  satisfied,  equation  (12)  may  be  solved  in  the 
^  with  the  functions 


neighborhood  of  P 

|u1(»1),©1(<i>1),u2(<p1),e2(»1),«p°(<i'1)|  e  C1 


(15) 


The  function  represents  the  ideal  law  of  motion.  In  most  cases 

(for  conjugate  tooth  action)  function  2 ^  1  ^  ^near* 

Equations  (6)  and  (7)  also  yield  a  system  of  five  independent  equations  in 
six  unknowns  (upOj^pU^,^,^) 

9  j  ( »9  ^ » ©2  »®2  2  ’  ^  =  0  (j  =  1»2,...,5)  (16) 

It  is  assumed  that  this  system  is  satisfied  by  a  set  of  parameters 

p(2>  .  (u<2>,«<2>,»<1>,42>,42>,42>)  a?) 


with  the  same  value  of  as  in  the  set  P^, 

P^  the  Jacobian 


If  in  the  neighborhood  of 


^  (9]^  >9£  *93  *94  »9g) 


*  0 


D(Up«p  u^  »©2  2  ^ 

then  system  (16)  may  be  solved  with  the  functions 

|  u1(»1AQ)se1(<i>1,AQ),u2(91,AQ),e2(<p1,AQ),<i>2(»1,aQ)|  e  C1 


(18) 


(19) 
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Function  ^(^.aQ)  represents  the  actual  law  of  motion  transformation  - 


the  law  of  transformation  of  motion  which  corresponds  to  errors  of  manufac¬ 
turing  and  assembly.  Kinematic  errors  of  the  gear  drive  are  represented  by 
the  function 

2 (•>> i » AQ )  *  v 2 ^  1 )  (20) 

This  method  of  solution  can  provide,  not  only  the  kinematic  errors  of  a  gear- 
set,  but  also  the  new  path  of  the  contact  point  (see  functions  (10)). 

In  general,  the  numerical  solution  of  a  system  of  five  nonlinear  equations 
is  a  difficult  problem  which  requires  many  iterations.  To  save  computer  time 
an  effective  method  of  solution  was  recently  proposed  by  Litvin  and  Gutman 
[2],  The  principle  of  this  method  follows: 

The  system  of  equation  (16)  may  be  represented  as 


^ * ^2 »®2 2 1 ’^2 ’ ^  =  ^ 

(21) 

f  2  ( i  >  ®i  i » ^2  *  ®2 2  ’  ^  ^  1 ’^2  ’  ^  =  ^ 

(22) 

^3^ul,el,<l’l,u2,e2’<p2,^’^l,^2’A^^  =  ^ 

(23) 

^4(ul»epu2’e2,<,>2’A^^  =  ^ 

(24) 

^5(ul’el’u2’€2,(,>2,A^  =  ^ 

(25) 

Equations  (21)  to  (23)  are  determined  from  vector  equation  (6),  and 
equations  (24)  to  (26)  from  vector  equation  (7).  Here,  A  represents  the 
shortest  distance  between  the  axes  of  rotation  of  the  two  gears  and  and 
represent  the  axial  settings  of  the  gears  (Fig.  2). 

9 


Systems  S^(X^,Y^,Zj)  and  S^^.Y 2,Z2)  sfl0wn  in  fi9ure  2  are  rigidly 
connected  to  the  driving  and  driven  gears,  respectively.  Now  suppose  that  some 
points  M^(uj,Oj)  and  on  surfaces  and  are  chosen. 

With  a  set  of  known  parameters  (u^.e^, 02,62) ’  eclljat^ons  (?4)  and  (25)  become 
a  system  of  two  equations  in  two  unknowns  which  may  be  expressed  as 


F  t  (v  ^  2  ^  =  0 

(26) 

F2^r?2^  =  0 

(27) 

Upon  solving  for  ^  and  <p ^ *  one  checks  that 

satisfied: 

the  following  equations  are 

A  -  ( u^  ,©j^  ,<p  j , U2  ,©2 ,9 2  >  )  =  0 

(28) 

8  2  -  K2  ( u^  ,©2  ,*P  ^  *  u2 ’^2 2  *  ^  =  ^ 

(29) 

H2  -  K3(u1,©2,<p1,u2,©2,<('2»a())  =  0 

(30) 

where  A,  Hj,  Hg,  and  &Q  are  given  values. 

In  general,  the  solution  of  the  above  two  systems  of  equations  ((26)  and 

(27)  and  (28)  and  (30))  requires  an  iterative  procedure.  In  practice,  one  of 

the  four  variable  parameters  fixed,  and  the  other 

three  are  changed  such  that  the  two  equation  systems  are  satisfied. 

The  advantage  of  the  above  method  lies  in  the  ability  to  divide  the  system 
of  five  equations  ((25)  to  (30))  into  two  subsystems  of  two  and  three 
equations,  and  to  solve  them  separately. 

This  method  was  applied  to  investigate  the  sensitivity  of  Wildhaber- 
Novikov  gears  to  errors  of  center  distance  mounting  [1].  These  gears  are 
generated  by  two  rack  cutters  which  have  normal  sections  as  shown  in  Fig.  3. 


Surfaces  of  the  rack  cutters  are  in  tangency  along  a  straight  line  M  -  M‘ 


which  is  parallel  to  axis  Z,  and  passes  through  point  M  of  the  normal 

a 

section.  In  the  normal  section,  the  shape  of  each  rack  cutter  is  a  circular 
arc  of  radius  (i  =  I,  II).  The  location  of  point  M  is  defined  by  the 
parameter  e  =  eQ  (Fig.  3).  The  line  M  -  M*  generates  a  helix  on  the  gear 
tooth  surface  which  is  the  path  of  contact  points.  Although  the  procedure 
described  above  is  primarily  for  numerical  solution,  in  this  case  analytical 
results  were  obtained. 

The  investigation  showed  that  the  error  of  center  distance  aA  resultea 
in  the  change  of  location  of  the  contact  point  path.  The  new  location  of  the 
path  is  represented  by  the  equation 


sin  e 


n 

ii 


^  (pj  sin  e  +  r^)^  +  (pj  cos  e  sin  x)2  (31) 

Where  R^  is  the  radius  of  gear  1  cylinder  on  which  lies  the  helix  of  the 
new  contact  point  path.  Parameters  e,  pj ,  pjj,  and  bj  j  are  shown  in 
Fig.  3;  r^  is  the  pitch  cylinder  radius. 

With  aA  =  0, 


sin  e 


0  = 


(32) 


where  eQ  is  the  parameter  corresponding  to  the  desired  location  of  the 
contact  point  path. 
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KINEMATIC  RELATIONS  BETWEEN  PARAMETERS  OF  CONTACT  POINT  MOTIONS 
The  following  relations  are  the  basis  of  the  second  method  for  the 
determination  of  kinematic  errors  in  gear  trains.  As  stated  above,  the 
tangency  of  gear  tooth  surfaces  is  represented  by  equations  (4)  and  (5). 
Because  of  the  continuity  of  tangency  of  these  surfaces,  it  is  required  that 


i  i )  = 


( ^2  2  ^ 


(33) 


♦  /  1  \ 

t  v ui ,<©i  »«p i , )  =  n^t;(u2,«2,*2) 

Here  (i  =  1,2)  is  the  velocity  of  the  surface  contact  point  in  abso¬ 
lute  motion  (with  respect  to  the  frame);  is  the  velocity  of  the  tip  of 

the  surface  unit  normal  in  absolute  motion  (also  with  respect  to  the  frame). 
Henceforth,  rj^  is  designated  as  v^  and  n|^  as  n^. 

The  velocity  of  absolute  motion  may  be  represented  as  the  sum  of  two 
components:  (a)  the  velocity  of  transfer  motion  (together  with  the  surface) 
and  (b)  the  velocity  of  relative  motion  (with  respect  to  the  surface). 
Consequently, 


;<z)( 


(34) 


v*1*  - 
_abs 

♦  i'11. 

v(2)  - 
_abs 

v<?>  ♦  v<2> 

_tr  _r 

(35) 

;u) . 

_abs  " 

n(2)  - 
IJabs 

n<2>  .  ;<2> 

_tr  _r 

(36) 

For  a  surface  represented  by  a  vector-function 
1^(u.,ei,»i)  (i  =  1,2) 


(37) 
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and  a  surface  unit  normal 


n^(u. 

(i  =  1, 

2) 

lowing  comes  from  def 

initions 

(a)  and  (b) 

above 

vo>  _ 

ar^ 

*i 

* 

v(i)  _ 

ar<f> 

dui 
__ _  + 

arO 

dei 

!tr  - 

3*i 

it- 

_r 

aui 

ZT 

3ei 

it- 

an*1* 

. 

an((> 

du- 

an^ 

dei 

^  - 

- 

7 

;(i)  _ 

- 

1 

_ _  + 

- 

1 

_tr  " 

3<P  ^ 

W' 

_r 

3Ui 

ir 

30i 

dt 

(38) 


(39) 


(40) 


Transfer  velocity  may  also  be  determined  in  a  kinematical  way  by  supposing 

that  a  gear  with  surface  z ^  rotates  about  an  axis  that  does  not  pass 

through  the  origin  0.  of  coordinate  system  S^,  which  is  rigidly  connected 

to  the  frame.  Vector  ^  ^  is  the  vector  of  angular  velocity  of  the  gear's 
rotation.  Then, 


(*>.  JUxr^  +  R 


tr 


x  w' 


(41) 


where  r^  is  the  position  vector  drawn  from  origin  0.  to  the  contact 
point  on  the  tooth  surface  and  is  a  vector  drawn  from  0^  to  an  ar¬ 

bitrary  point  on  the  gear's  axis  of  rotation. 

The  transfer  velocity  is  represented  by  the  equation 


x  n 


(i) 


(42) 


Equations  (33)  to  (36)  yield  the  following  kinematic  relations  for  two 
tooth  surfaces  which  are  in  continuous  tangency: 


(43) 
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(44) 


nlD 

r 


- 


n<2> 

_r 


Equations  (43)  and  (44)  were  first  proposed  by  Litvin  [1],  On  the  basis  of 
these  equations,  important  problems  in  the  theory  of  gearing  were  solvea,  such 
as  avoiding  tooth  undercutting,  deriving  the  relations  between  curvatures  of 
two  gear  tooth  surfaces  in  mesh,  and  determining  the  kinematic  errors  of  gear 
drives  which  are  caused  by  errors  of  manufacturing  and  assembly. 

APPROXIMATE  METHOD  OF  CALCULATION  OF  GEAR  DRIVE  KINEMATIC  ERRORS 
As  a  general  rule,  kinematic  errors  of  a  gear  drive  determined  by  the 
exact  method  must  be  obtained  in  a  numerical  way  using  a  computer.  This  is  a 
disadvantage  of  the  exact  method.  Therefore,  an  approximate  method  with  the 
opportunity  to  obtain  accurate  results  analytically  is  now  presented. 

Figure  4  shows  two  gear  surfaces  z^  and  l2  which  are  not  in  tangen- 
cy  due  to  errors  of  manufacturing  and  assembly.  Points  and  do  not 

coincide,  position  vectors  and  rj2^  are  not  equal,  and  surface  unit 

normal  vectors  ni^  and  n|2  ^  do  not  coincide.  To  bring  the  two  surfaces 
into  contact  it  is  sufficient  to  hold  one  gear  fixed  and  rotate  the  other  gear 
by  an  additional  small  angle.  Since  the  gear  with  surface  z-^  is  the  ariv- 
ing  gear,  it  is  preferable  to  fix  the  position  of  surface  z^  and  r°tate 
surface  A<p ,,  to  bring  it  back  into  contact  with  z^.  The  additional 
angle  of  rotation  Aq>2  represents  the  change  of  the  theoretical  angle  of 

rotation  <p°  which  is  exerted  by  errors  of  the  manufacturing  of  manufacturing 
and  assembly.  The  A<p ^  is  as  yet  an  unknown  function  of  the  vector  of  errors 
aQ  and  varies  in  the  process  of  motion.  Thus 

A i*2  =  ftfj.AQ)  (45) 


1 .  jiaSAjSssfeV  ■ 
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The  determination  of  the  Function  (45)  is  based  on  the  following  kinematic 


relations,  which  are  analogous  to  (43)  and  (44): 


ds<‘)  *  ds<‘>  *  dsO)  .  ds<2>  ♦  ds<2>.  ds<2> 
_tr  _r  __q  _tr  _r  _q 


dnj1)*  dn<H  *  dn‘l>  .  dn<2>  .  d„t2>  *  dn<2> 
„tr  „r  _q  _tr  _r  _q 


where  ds^  is  the  displacement  of  contact  point  and  dn^  (i  =  1,2)  is 
the  change  in  direction  of  the  surface  unit  normal  due  to  errors  of  manufac¬ 
turing  and  assembly.  To  bring  the  surfaces  into  contact,  it  is  sufficient 

to  rotate  only  gear  2,  holding  gear  1  at  rest.  Therefore,  ds|^,  and  anj^ 
are  zero,  and 

ds)1*  ♦  ds*1’  .  ds<2!  ♦  ds<2>  *  ds<2>  (4t 


dn^1^  +  dn^  =  dn|2^  +  dnj.2^  +  dn^ 
_r  _q  _tr  _r  „q 


To  determine  relations  between  dsi^,  ds^,  and  dsj^  take  the  follow¬ 
er  __q  _q 

ing  scalar  products: 

n  .  (dsi1)  +  ds^)  =  n  .  (ds^5  +  ds^2)  +  dsj,2))  (5C 


Since  vectors  dsj^  and  dsj.2^  must  lie  the  common  tangent  plane  T,  equation 
(50)  is  reduced  to  the  following: 

n  .  ds[2^  *  n  .  (ds^  -  ds^)  (51) 


The  vector  ds^'  may  be  represented  by  the  following  cross  product: 
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(52) 


where  d?^  is  the  incremental  angle  of  rotation  of  gear  2  and  r^  is  the 
position  vector  drawn  from  an  arbitrary  point  of  the  axis  of  rotation  to  the 
contact  point  M. 

Equations  (51)  and  (52)  yield 


(53) 


Equation  (53)  is  the  basic  equation  for  the  determination  of  kinematic  errors 
of  gear  drives.  Its  application  will  be  demonstrated  in  the  following  sections. 

Analogous  scalar  products  can  be  composed  on  the  basis  of  equation  (49). 

It  can  be  proven  that  these  scalar  products  are  zero  because  the  vectors  in 
equation  (49)  all  belong  to  the  tangent  plane.  Hereinafter,  the  following 
notations  will  apply: 


(54) 


where  EAq^  and  z Aq^represent  the  sum  of  linear-error  vectors  due  to 
manufacture  and  assembly  of  gears  1  and  2,  respecti vely. 

In  many  cases,  however,  errors  in  gear  trains  do  not  result  from  linear 
displacements,  but  rather  from  angular  displacements.  For  instance,  kine¬ 
matic  errors  may  result  from  the  misalignment  of  gear  shafts. 

Figure  5  shows  the  axis  of  gear  2  rotation  a-a  in  its  ideal  position. 
Suppose  that,  due  to  an  error  of  assembly,  axis  a-a  is  rotated  about  a 
nonintersecting  axis  B-B.  Such  an  error  of  assembly  may  be  represented  by 
the  vector  &«,  which  is  directed  along  axis  B-B,  where  the  direction  of 
corresponds  to  the  direction  of  rotation  by  the  right-hand  rule. 
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With  the  given  vector  ts  the  displacement  of  contact  point  M  may 

be  determined  as  follows: 

(a)  Vector  aa,  directed  along  the  axis  B-B,  is  replaced  by  an  equal  vector 
aa,  which  passes  through  the  origin  0^  and  the  vector-moment  R  x  a«.  Here 

R  is  a  position  vector  drawn  from  0^  to  an  arbitrary  point  on  the  line  of 
action  of  vector  aa  (Fig.  5). 

(b)  The  displacement  aqv  '  corresponding  to  aa  may  be  represented  by 

aq^  =  aa  x  r./M^  +  R  x  aa  =  aa  x  (r£M^  -  R)  (55) 

A  similar  equation  may  be  developed  to  determine  the  displacement  of  the  contact 
point  M  exerted  by  an  angular  erro1'  corresponding  to  gear  1. 

With  notations  (54)  the  equation  (53)  for  the  determination  of  kinematic 
errors  may  be  represented  as  follows: 

(a x  r^  +  zaq)  .  n^  =  0  (56) 

where  zaq  =  zaqi.2^  -  zaq^  and  n^M'  is  the  surface  unit  normal  at  the  contact 
point  M. 

The  location  of  the  contact  point  M  and  the  direction  of  the  unit  normal 
n^  change  in  the  process  of  motion.  A  further  simplification  of  equation 
(56)  results  by  assuming  that  in  all  positions  the  contacting  tooth  surfaces 
have  a  common  normal  which  passes  through  the  contact  point  M  and  the  pitch 
point  P.  This  is  the  fundamental  law  for  uniform  motion  transmission.  For 
planar  gears  the  pitch  point  P  coincides  with  the  point  of  tangency  of  the 
pitch  circles  (gear  centrodes).  The  pitch  point  for  bevel  gears  is  located  on 
the  line  of  tangency  of  pitch  cones.  In  both  cases  the  surface  unit  normal  n 


is  col  linear  to  PM  and 


r<">  r<P>  *  W(  ,  r<P>  *  *„<"> 
Equations  (56)  and  (57)  yield 

(*(2>  xr<P>  +£Aq)  .  nW-0 


(57) 


(58) 


because 


U*<2>  x  r<">)  .  n<">  . 


A</2^  x  (riP^  +  xn^) 


-  (a*<2>  x  r<P>)  .  nW  + 


n<M>) 


(2)  (M)  (M) , 

Aq> '  'xnv  'nv  ' ) 


U./2)  X  r<P>  )  .  nW) 


(59) 


Application  of  equation  (58)  in  place  of  equation  (56)  has  the  advantage 

that  the  location  of  the  pitch  point  may  be  considered  as  a  constant  (r^ 

=  const).  However,  the  direction  of  the  surface  unit  normal  is  a  function  of 
2 •  Three  types  of  gears  -  involute  (spur  and  helical)  and  Wildhaber- 
Novikov  -  are  exceptions  to  this  statement.  For  these  gears  the  unit  normal 
of  the  gear  surfaces  at  their  contact  point  does  not  change  its  direction. 

8ecause  of  kinematic  errors,  the  angular  velocity  ratio  fluctuates  as  the 
gear  teeth  pass  through  mesh.  Figure  6  shows  functions  for  two  types  of  kine¬ 
matic  errors.  The  first  is  a  piecewise,  nonlinear,  periodic  function  which 
has  a  period  that  depends  on  the  ratio 


"12 


"7 


(60) 


where  (i  =  1,2)  is  the  number  of  gear  teeth  and  b  and  a  are  the 
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minimum  integral  numbers  with  which  the  ratio  can  t>e  expressed.  The 
angle  of  rotation  of  gear  1  corresponding  to  the  period  of  function  A<p'  '[ <p 
is  equal  to  2*a. 

Such  functions  of  kinematic  errors  are  caused  by  (a)  the  eccentricity  of  gears  ana 
(b)  the  crossing  of  the  theoretical  axis  with  the  axis  of  rotation  of  gears.  (The 

shortest  distance  between  these  axes  rotates  in  the  process  of  motion.) 

(21 

The  second  type  of  function  Abv  '(b,),  shown  in  Fig.  6(b),  has  a  period  of 

2tt  ^ 

of  bj  =  This  function  is  exerted  by  (a)  errors  in  the  generating  process 

of  gear  teeth  and  (b)  errors  of  gear  axis  location  which  do  not  change  in  the 
process  of  meshing,  etc. 

APPLICATION  OF  THEORY  TO  ECCENTRICITY  OF  INVOLUTE  SPUR  GEARS 

Figure  7  shows  base  circles  of  radii  r^  and  rj^  for  two  involute 
spur  gears.  The  rotation  centers  of  the  gears  are  denoted  0^  and  02* 

If  the  centers  of  base  circles  0^  coincide  with  centers  of  roation  0^ 

(i  =  1,2),  and  then  vectors  of  gear  eccentricity  Ae^  =  O^oT  are  zero.  The 
involute  curves  are  in  tangency  at  a  point  M  of  the  line  of  action  KL. 

To  model  the  meshing  of  gears  with  eccentricity,  gears  1  and  2  are  trans¬ 
lated  from  their  theoretical  positions  by  =  Ae^  (i  =  1,2).  Now,  the 

center  0^  will  be  offset  from  the  center  of  rotation  0^'.  Because  of 
this  displacement  of  the  gears,  the  tangency  of  their  involute  curves  is  brok¬ 
en:  the  curves  will  wither  interfere  with  each  other  (intersect)  or  lose  con¬ 
tact.  To  bring  the  involute  curves  into  the  contact  once  again,  it  is  suffic¬ 
ient  to  rotate  gear  2  by  a  small  angle  a</^.  According  to  equation  (68) 

(2) 

the  angle  Atp  may  be  determined  with  the  equation 


(62) 


where  b0  is  the  pressure  angle. 

Vectors  of  eccentricity  ae^  and  ae2  form  angles  and  b2  with  vector 

0(1>0<2>;  these  angles  are  measured  in  the  direction  of  gear  roatation  (Fig.  '/). 
The  dot  products  yield 

,(M) 


ae,  .  n' 


=  aej  sin(Bj  +  *Q) 


,(M) 


ae2  .  nv  '  =  ae0  sin(+0  -  b2) 

It  results  from  equations  (61)  to  (62)  that 

(2)  aex  sinUj  +  ♦Q)  +  ae2  sin(s  -  *0) 


abv 


m 


(63) 


(64) 


where 


r ^  =  r ^  cos  ii»0 


(65) 


is  the  radius  of  the  base  circle  of  gear  2. 

The  center  (i  =  1,2)  of  the  base  circle  rotates  in  the  process  of  meshing; 
oj1*  and  oj^are  two  instantaneous  positions  of  this  center  (fig.  8).  Angles 
and  &2  can  be  represented  as  follows: 

'r 


ei  =  B10  +  »■ 


b2  “  e20  +  ’2 


(66) 


where  6^  and  b2q  correspond  to  the  initial  positions  of  centers  0^  and  0 2> 
with  =  (j>2  =  0. 


Equations  (64)  and  (66)  yield 

(2)  aej  sin(?1  +  v±)  +  ae2  sin(<p2  +  y2) 


A<f> 


.(2) 


(67) 


where 

Y1  =  {810  +  *o)  ;  y2  =  {e20  -  *o] 

For  convenience,  consider  the  kinematic  error  function  to  have  zero  magnitude  at 
=  0.  Then,  the  kinematic  error  becomes 
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(2)  =  w(2)(n)  -  av(2)(0) 

Ae^  [sin(b1  +  y^)  -  sin  y^I 


xrr 


»n  + 


ae2  [sin(b?  +  y2)  -  sin  ^ 


JTT 


(68) 


where 


"21 


JD  M 

^2  _  ^b _ ^1 

~  r(2)  "  N 

“1  rb  N2 


Equation  (68)  represents  the  kinematic  error  of  a  gear  train  with  two  gears  as 
the  sum  cf  two  harmonics.  The  periods  of  these  harmonics  are  equal  to  the 
periods  of  complete  revolutions  of  the  gears. 

Equation  (68)  may  be  made  symmetric  as  follows: 


(2)  Ae-, [ s i n(<p +  y^ )  -  sin  yj]  Ae2C s i n(«p 2  +  y2)  ~  sin 

~  — ± - 1 - i - m01  + - - - - 

jo  21  d2> 


"22 


2  Ae,[sin(b1  +  y^)  -  sin  Yi3 
i=l 

rb 


m2i 


(69) 


Here,  m22  =  1  and  <P 2  =  ’l  m21 

Equation  (69)  can  be  generalized  for  a  train  with  n  gears  as  follows: 


J  aei[sin(?.i  +  y,)  -  sin  y1- 3 

i-i - rrn 

rb 


m  . 


ni 


(70) 


where  a«^  is  the  resulting  kinematic  error  of  the  gear  train  represented 

as  the  angle  of  rotation  of  gear  n  (the  output  gear). 

A  complicated  gear  train  is  a  combination  of  pairs  of  gears.  The  parameter 
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y^  may  be  represented  as 


:  t 


for  the  driven  gear  of  the  pair.  For  instance,  for  computational  purposes,  a 
train  of  three  gears  must  be  replaced  by  two  pairs  of  gears.  The  idler  (inter¬ 
mediate  gear)  is  considered  as  the  driven  gear  in  the  first  pair,  and  as  the 
driving  gear  in  the  second  pair. 

Designate  the  kinematic  error  exerted  by  the  eccentricity  of  gear 

number  i  as 


Aei  " 


Aei[sin(<*>  ■  +  Yi)  -  sin  Yil 

Ji) 


(71) 


where  is  the  error  of  the  rotation  angle  9,. 

possible  value  of  this  error  is 


The  maximum 


AO-  -AO- 

i,max  i.min 


2  Aei 

loT 


(72) 


The  kinematic  error  of  the  train  may  be  represented  as 


ao^  =  z  AO-  m  . 

.  ,  i  ni 
i=l 


(73) 


Usually  gear  trains  are  applied  for  the  reduction  of  angular  velocities 
and  thus  mpi  is  less  than  1.  It  results  from  equation  (73)  that  the  last 
gears  of  a  train  (numbers  n,  n  -  1,  n  -  2)  induce  the  largest  part  of  the 
resulting  kinematic  error  ao^.  Therefore,  the  precision  of  these  gears 
must  be  higher  than  the  others. 
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The  largest  value  of  the  kinematic  error  function  ae^  and  its  distri¬ 
bution  above  and  below  the  abscissa  depend  on  the  combination  of  parameters 
y j  (i  =  1,  2,  ....  n).  figure  9  shows  the  distribution  of  a  function  ae^ 

(v.j)  exerted  by  eccentricity  of  gear  i  of  the  train. 

The  resulting  errors  of  a  gear  train  may  be  compensated  for  in  part,  by 
definite  rules  of  assembly  of  gears  with  eccentricity.  For  instance,  for  gears 
with  tooth  numbers  N^  =  N,>  and  equal  eccentricities  ae^  =  ae2  the 
resulting  kinematic  error  will  be  approximately  zero  if  eccentricity  vectors 
ae^  and  ae^  (fig.  7)  are  directed  opposite  each  other. 

APPLICATION  Of  THEORY  TO  ECCENTRICITY  OF  SPIRAL  BEVEL  GEARS 
For  spatial  gears  the  word  "eccentricity"  is  used  to  describe  that  the 
geometric  axis  of  a  gear  is  parallel  to,  but  does  not  coincide  with,  its  axis 
of  rotation  (Fig.  10).  As  the  eccentric  gear  rotates  its  geometric  axis  gener¬ 
ates  a  cylindrical  surface  of  radius  ae.  The  eccentricity  vector  ae  is  a 
vector  which  rotates  about  the  gear  axis.  The  initial  position  of  vector 
ae  (its  position  at  the  beginning  of  motion)  is  given  by  angle  a  (Fig.  10). 

Figure  11  shows  coordinate  systems  S^(XpYpZ^)  and  Sf(Xf ,Yf ,Zf ), 
which  are  rigidly  connected  to  gear  1  and  the  frame,  respectively.  System 
is  an  auxiliary  coordinate  system,  which  is  also  rigidly  connected  to 
the  frame.  Driving  gear  1  rotates  about  axis  Z^.  The  position  of  ae^ 
in  coordinate  system  is  given  by  the  angle  aj,  which  is  made  by 
aej  and  axis  X^.  The  current  position  of  ae^  in  coordinate  system  Sf 
(or  S^)  is  defined  by  the  angle  (<p^  +  a^)  and  the  matrix  equality 

Ue)1’]  -  [lfh] 


cos  y ^  0  sin  y j 

ae1  cos^  +  o^) 

0  1  0 

-ae^  sin((p^  +  o^) 

-sin  yj  0  cos  yj 

- 1 

o 

_ 1 

tV*  » 
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Here  [Ae^]  is  the  matrix  of  vector  Ae^  in  terms  of  its  projections  on 
on  the  axes  of  coordinate  system  S^.  The  3  by  3  matrix  [L^]  trans¬ 
forms  elements  of  the  column  matrix  Ae^  _to_  coordinate  system 
from  coordinate  system  Sh- 

Matrix  equality  (74)  yields 


]  = 


Ae^COS(<p^  +  ) COS  YJ 

-Ae1  sin(?1  +  a^) 
-Ae^cos((f>^  +  a^)sin  y^ 


(75) 


The  vector  of  eccentricity  of  the  driven  gear  can  be  defined  in  a  similar 
way.  Figure  12  shows  coordinate  systems  and  rigidly  connected  to 
gear  2  and  the  frame,  respectively.  The  auxilliary  coordinate  system  is 
also  rigioly  connected  to  the  frame. 


Vector  Ae^  is  represented  by  matrix  equalities 

u42>]  -  [Lfp] 


cos  y2  0  -sin  y2 

ae2  cos(<p2  +  (*2) 

= 

0  1  0 

Ae2  sin((j>2  +  a2) 

sin  y2  0  cos  y2 

0 

(76) 


which  after  matrix  multiplication  gives 


r  Ae2  cos(<j>2  +  a 2)  cos  y2 

Ue<2)]  = 

ae2  sin(<?2  +  a^) 

Ae2  cos(92  +  a2)  sin  y2 

(77) 


Kinematic  errors  induced  by  gear  eccentricities  may  now  be  found  by  applying 
Equation  (58)  as  follows: 


*t2M2V2) 


ae^  -  ae^ 


L 


,(2) 


(78) 


Here  aej^  and  ae^  are  given  by  matrices  (75)  and  (77);  vector  a<p|2' 
(Fig.  12)  is  represented  by  the  matrix 

l»<2>:  -  [Lfp]  t«<2>] 

— %  r— 

0 
0 

a<p 2 


cos  y2  0  -sin  y2 

0  1 

0 

sin  y2  0  cos  y2 

L_ 

-&<p  2  s  i  n  y  2> 

0 

aip  2  cos  y2 

__J 

(79) 


(2) 

Vector  '  represents  the  position  vector  of  the  contact  point  which 
belongs  to  the  line  of  action  and  n^  represents  the  common  unit  normal  of  the 
gear  surfaces  at  their  point  of  tangency. 

Equations  (78)  and  (79)  yield 


n  zae  +  n  zae  ♦  n  rae 

a*2  =  - * - 5 - 1 - 1 - i - i -  (80) 

-y  cos  y2  +  (x  cos  +  z  sin  Y2)n^  -  y  sin  y2  nz 

where  zae  =  ae^  -  aef2^,  zae  =  ae^  -  ae^,  zae  =  ae^  -  ae^  and 
xx  x  y  y  y  zz  z 

the  subscript  f  was  dropped. 

Projections  n  and  n  of  the  surface  unit  normal  n,  and  coordinates 
y  z 

x,  y,  and  z  of  the  contact  point  change  in  the  process  of  motion.  But  since  the 
changes  in  these  variables  are  relatively  small,  they  may  be  neglected  (Fig.  13). 


nf  =  sin  cc  if  +  cos  c(  (cos  6  jf  +  sin  6  kf) 
*f  =  0,  yf  =  0,  zf  =  L 


(81) 
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Equations  (80)  and  (81)  yield 


*2(*,)  -  A(*l} _ 

L  sin  *2  cos  *  cos  B 


where 

A(<p^)  =  a1  sin  (^  +  cos  (<?1  +  o^)  + 

&2  sin  (<?2  +  a 2 )  +  b2  cos  (^  +  <*) 
a-^  =  -ae^  cos  9c  cos  6 
a2  =  -ae2  cos  4>c  cos  b 

b1  =  ae^cos  y^  sin  4>c  -  sin  y^  cos  4»  sin  b) 
b2  =  -ae2(cos  y2  sin  *c  +  sin  y2  cos  *c  sin  b) 

N1 

v2  "  *1  N2 


(82  J 


It  is  concluded  from  the  form  of  equation  (82)  that  kinematic  errors  in¬ 
duced  by  the  eccentricity  of  spiral-bevel  gears  may  be  represented  as  the  sum 
of  four  harmonics:  the  period  of  two  harmonics  coincides  with  the  period  of 
revolution  of  gear  1;  the  period  of  the  other  two  coincides  with  the  period  of 
revolution  of  gear  2. 

The  function  a*2(<p.)  represented  by  equation  (82)  is  a  smoothed, 
continuous  function  which  serves  as  a  first  approximation.  In  reality  the  true 
function  ^ j_ )  breaks  as  different  sets  of  teeth  come  into  mesh.  This 
break  can  be  discovered  if  a92(<P^)  is  determined  by  equation  (80). 

GEAR  TRAIN  VIBRATION  AND  NOISE  MEASUREMENT 
To  illustrate  the  principles  discussed  on  the  subject  of  gear  train  preci¬ 
sion,  Figs.  14  and  15  are  used.  These  figures  show  some  frequency  spectrum 
measurements  made  on  a  helicopter  transmission  running  in  a  test  stand  [7], 

The  transmission  had  a  spiral-bevel  input  stage  with  19  teeth  on  the  pinion 
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and  71  teeth  on  the  gear.  The  pinion  was  turning  at  6200  rpin  and  the  output 
shaft  at  365.5  rpm.  The  output  stage  was  a  spur  planetary  arrangement  with  a 
27-tooth  sun,  3  planet  gears,  each  with  35  teeth  and  a  yy-toolh  ring  gear  which 
was  splined  to  the  transmission  housing.  An  accelerometer  was  mounted  on  the 
case  immediately  outside  the  spline. 

Figure  14  shows  a  broadband  frequency  spectrum  measurement  of  the  vibra¬ 
tion  signal.  The  spur  mesh  frequency  was  583  Hz  and  the  spiral  bevel  mesh 
frequency  was  1963  Hz.  The  spiral-bevel  vibration  signature  was  much  stronger 
than  the  spur  signature.  This  indicates  that  the  meshing  accuracy  was  better 
for  the  spur  mesh  than  for  the  spiral-bevel  mesh.  There  are  also  other  peaks 
in  the  spectrum  at  multiples  of  the  fundamental  frequencies  of  1963  Hz  and  583 
Hz.  These  other  peaks  are  the  higher  harmonics  due  to  the  noise  and  vibration 
pulsations  as  the  teeth  mesh  being  different  from  the  pure  sinusoidal  shape  as 
shown  in  Fig.  9. 

Figure  15,  an  expanded  region  of  the  autospectrum  plot  given  in  Fig.  14, 
shows  many  peaks  that  are  symmetrically  located  about  the  spur  gear  mesh  fun¬ 
damental  frequency  peak  at  583  Hz.  These  peaks  locate  the  sideband  frequen¬ 
cies  which  are  due  to  sources  of  modulation  in  the  time  dependent  vibration 
waveform.  Each  source  of  modulation  may  produce  one  pair  of  sidebands  if  it 
is  a  harmonic  modulator.  If  nonharmonic,  the  side  bands  will  repeat  many 
times,  as  in  the  case  in  Fig.  15. 

SUMMARY  OF  RESULTS 

Kinematic  precision  is  affected  by  errors  that  are  the  result  of  either 
intentional  adjustments  or  accidental  defects  in  the  manufacturing  and  assem¬ 
bly  of  gear  trains.  A  general  method  for  the  determination  of  kinematic  pre¬ 
cision  of  gear  trains  has  been  explained.  The  general  method  is  based  on  the 
exact  kinematic  relations  for  the  contact  point  motions  of  the  gear  tooth  sur- 
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faces  under  the  influence  of  errors.  An  approximate  method  was  also  explained. 
Example  applications  of  the  general  and  approximate  methods  were  demonstra¬ 


ted  for  gear  trains  consisting  of  involute  (spur  and  helical)  gears,  circular 
arc  (Wi Idhaber-Novikov)  gears,  and  spiral-bevel  gears.  Gear  noise  measurements 
from  a  helicopter  transmission  were  presented  and  discussed  with  relation  to 
the  kinematic  precision  theory.  The  following  results  were  obtained: 

1.  The  exact  numerical  iterative  procedure  for  finding  kinematic  errors, 

Ab£,  is  as  follows:  From  equation  (20)  find 

a< p2  =92(9^  aQ)  -  92(9^ 

where  the  angles  9 ^  ana  ^  have  ^een  determined  from  an  iterative 
solution  of  the  nonlinear  algebraic  system  of  five  equations,  which  are  separ¬ 
able  into  two  systems  of  two  and  three  equations  each  as  follows: 

i" } (9 1 »  92)  =  0,  i  =  1,2 
^  A  -  Kj(Uj,Op9pU2,©2>92»AQ)  =  0 

^  ^1  -  l<2 (  U^,8p9 -^2,02*9 2 »AQ)  =  0 

^  ~  K|j(  Oj  ,6j  ,9  j  ,  U2  ,©2 ’9g  >  aQ  )  =  0 

where  A,  Hj,  Hg,  and  aQ  are  given  values  and  9^  is  determined  by  solving  with 
aQ  =  0,  whereas  92  is  determined  by  the  full  solution  of  the  five  equations. 

2.  The  approximate  equation  for  kinematic  error  A9^  is 

Ubl2>  x  r<P>  x  S.,)  •  n<M>  .  0 

where  r ^  is  the  radius  vector  to  the  pitch  point,  n^  is  the  surface  normal 
vector  at  the  contact  point  M,  and  £aq  is  the  sum  of  known  error  vectors. 


3.  Application  of  the  formulas  showed  that  Wi Idhaber-Novikov  gears  are  sen¬ 
sitive  to  any  errors  which  cause  changes  in  the  center  distance.  A  formula  for 
the  location  of  the  gear-tooth  contact  point  was  given. 

4.  It  was  found  that  for  a  pair  of  spur  gears  the  kinematic  error  function 
due  to  eccentricities  is  a  sum  of  two  simple  harmonics.  For  a  multistage  speed 
reducer,  it  was  concluded  that  accuracy  in  the  final  stages  has  the  most  impact 
on  kinematic  error.  For  gears  with  approximately  equal  known  eccentricities, 
the  kinematic  error  may  be  compensated  for  by  directing  the  eccentricities 
opposite  one  another. 

5.  For  a  pair  of  spiral-bevel  gears,  the  kinematic  error  function  due  to 
eccentricities  is  a  sum  of  four  harmonics. 
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